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Abstract 

The possibility that the observed Ultra High Energy Cosmic Rays 
are generated by high energy neutrinos creating "Z-bursts" in reso- 
nant interactions with the background neutrinos has been proposed, 
but there are difficulties in generating enough events with reasonable 
incident neutrino fluxes. 

We point out that this difficulty is overcome if the background 
neutrinos have coalesced into "neutrino clouds" — a possibility pre- 
viously suggested by some of us in another context. The limitations 
that this mechanism for the generation of UHECRs places on the high 
energy neutrino flux, on the masses of the background neutrinos and 
the characteristics of the neutrino clouds are spelled out. 
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1 Introduction 



The continuation of the cosmic ray spectrum beyond the GZK (Greisen, 
Zatsepen and Kuzmin) cutoff is a puzzle which has inspired many attempts 
at explanation. Neutrino-anti-neutrino annihilation into a Z-boson which 
subsequently decays to hadrons and leptons is one promising solution of the 
problem. Should relic neutrinos possess a mass in the eV/c 2 range the center 
of momentum energy required for resonant annihilation can be produced by 
cosmic ray neutrinos of energy ~ 10 21 eV. This proposal was first examined in 
detail by Weiler [[J and Fargion @ , who assumed that the relic neutrinos 
cluster in galaxies up to two or three times the mean relic density, and found 
that the required flux of cosmic ray neutrinos is larger than had previously 
been suggested. It has been suggested that the clustering may extend to 
densities of 100 to 10000 times the normal relic density on galactic scales, to 
moderate the required incident flux ^j. 

Independently, it has been shown that the relic neutrinos may have a 
local density, on scales of the solar system to a few parsecs, which may be 
many orders of magnitude higher than standard cosmology predicts with- 
out conflicting with any known experimental results ||. Furthermore, the 
formation mechanism provides a strong argument for a natural association 
between the so called neutrino clouds of Ref. (0) with small celestial objects 
such as stars and solar systems. 

In this note we investigate the consequences of this latter scenario for the 
measured Z-burst rate and for the bounds on the neutrino cloud. Specifically, 
we find the neutrino cloud radius and density range required to reproduce the 
observed flux of UHECRs and use the non-observation of coincident UHECR 
events to strengthen further these constraints. 

We note that recent studies || [| of the Z-burst explanation for UHECRs 
have derived information about the masses of the relic neutrinos. In the light 
of the present study of a dense local cloud of neutrinos as the target, we note 
that the Fermi motion of the neutrinos in the local cloud would complicate 
attempts to determine the neutrino mass spectrum from the observations of 
UHECRs. 

2 Numbers and formalism 
2.1 Z-burst flux 

To evaluate the Z-burst flux at earth we will follow closely the formalism used 
by Weiler ||. The change in the differential neutrino flux density F u (E u ,x) 
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at a given energy E v and a distance x from its source is given by 

dF v (E v ,x) = -F v {E v ,x)a v {E v )n v {x)dx ) (1) 

where <j v (E v ) is the annihilation cross section, and n v (x) is the local neutrino 
density. The reduction in neutrino flux found at a distance D from the source 
is therefore found to be 



5F V (D) = J dE u F u (E u ,0)\l-exp{-<T(E u )S(D)}\, (2) 

where S(D) = f®n v (x)dx is the column density, and 5F U (D) = F u (0) — 
F V (D). 

In the narrow resonance approximation, where the variation in F U (E U , 0) 
over the resonant energy is small, we can write 

6F V (D) = E R F v {E a ,0)J^(l-ew{-a(a)S(D)}\ (3) 

where the resonance energy Er corresponds to a neutrino energy such that 
the center of momentum energy is s = M|. From this point on references 
to the neutrino flux should be taken to mean the quantity ErF u (Er, 0) with 
units cm _2 s~ 1 sr~ 1 . The main contribution to a(s) is given by the s-channel 
Z-resonance, vv — > f f, with a branching ratio of 70% to hadrons, 20% to 
neutrinos and 10% to charged leptonspf, 
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f 



m)-2t 3 (f)Q f s 2 w + 2Q}s* 



f w 



[s - M\Y + Y\M\ 



(4) 



where the quantities in the square brackets are the standard model weak 
isospin numbers, charge and sine of the Weak angle. In our calculation we 
set a(s) = cr z (s). 

The only quantity entering these considerations which has been measured 
directly are the post GZK events, with fluxes quoted = 10 19 cm 2 s 1 sr 1 

above 5 x 10 19 eV and F ph ~2x lO^cm-V^r- 1 above 5 x 10 20 eV @@. 
Assuming the Z-burst flux equals 5F U (D), and that the Z's are then converted 
to nucleons and photons, we can write 

F ph = M x 5F V (D), (5) 

where M is the photon and nucleon multiplicity per Z-burst, estimated to be 
~ 30. The aim of any calculation must be to reproduce the measured value 
of F p / 7 . 
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2.2 Neutrino flux 



The Z-burst model, be it in the presence of a non-clustered relic neutrino 
background or a neutrino cloud, requires a flux of UHE neutrinos. It is 
no surprise that the neutrino flux has not been directly measured at UHE 
energies; however there are theoretical as well as experimental upper bounds. 

The experimental upper bounds are derived in Ref. (0) from the non- 
observation of deeply penetrating particles with energy > 10 17 eV at the Fly's 
Eye detector in Utah [IU|. The upper bounds on the flux are: 2 x 10~ 13 , 3 x 
KT 14 and 5 x KT 16 cnrV^- 1 for E v = 10 17 , 10 18 and 10 20 eV respectively. 
This is unfortunately not able to constrain the Z-burst model. In the absence 
of experimental data one needs to turn to models of astrophysical particle 
acceleration. 

The mechanism by which primaries with energies above 10 19 eV are pro- 
duced is as much a puzzle as their propagation to the earth. Neutrino pro- 
duction resulting from the interaction of nucleons produced in AGN jets or 
gamma-ray bursts interacting with the cosmic microwave background is one 
explanation; other explanations involve more exotic physics such as decaying 
super heavy particles, or topological defects. For a review of such theories 



and a discussion of the expected neutrino flux see Ref. (fill) and references 
therein. The former explanation is the only theoretical description which is 
non-speculative in that it uses only standard model physics. There is however 
an ongoing debate as to the magnitude of the neutrino flux which may be 
derived from conventional sources [ 12] . Depending on the assumptions made 



the upper-bound to the neutrino flux from astrophysical sources ranges from 
E R F„(E R ,0) < 3 x KT 20 to 1 x 10- 19 cm- 2 s- 1 sr- 1 . 



2.3 Cloud density and radius 

While the neutrino cloud can be formed with a complex density profile, e.g. 
with more massive species at the center ||, we parameterise the cloud as a 
sphere with no diffuse boundary such that the density of the relic neutrinos 
is 



n v {x) = n c (x)9(R — x) + n (x)9(x — R), (6) 

where n c (x) is the density inside the cloud and n (x) is the density outside. 
In the standard theory the density of relic neutrinos is no(x) = 54cm -3 ; 
the formation of neutrino clouds would deplete this number to a negligible 
quantity, such that in our calculations it can be safely ignored. A more 
complex model of the cloud structure would include a sum over the different 
mass eigenstates of the neutrino species, as well as a diffuse boundary, and 
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perhaps some directional and temporal dependence taking into account cloud 
dynamics. However these additions will not result in large changes to the 
calculated flux. We consider relic and cosmic ray neutrinos with a mass of 
leV/c 2 , and clouds with radii in the range of 10 14 — 10 20 cm [] and density in 
the range of 10 10 — 10 16 cm -3 . These parameters are not in conflict with the 
requirement that the total mass of dark matter within the orbit of Uranus 
be less than ~ 3 x l(r 6 M pi . 



3 Results for neutrino clouds 

Clouds of neutrinos with density ranging from 10 10 cm" 3 to 10 16 cm -3 have a 
correspondingly large range in Fermi energy. In fact 1 eV/c 2 mass neutrinos 
in a 10 10 cm -3 cloud are non-relativistic, while in a 10 16 cm -3 cloud they are 
very much relativistic. In the following analysis we temporarily neglect the 
effects of the Fermi motion for clarity in the exposition, returning to it at a 
later point. 

The crucial element in our analysis is the column density S(D). In Fig. p] 
the column density needed to produce F p / 7 for various values of neutrino flux 
is shown. The upper curve corresponds to a flux of F p / 7 = 10~ 19 cm~ 2 s _1 sr _1 
while the lower curve corresponds to F p / y = 10~ 20 cm -2 s -1 sr~ :L . We would 
expect the actual value lies between these two values. The minimum col- 
umn density for a source of neutrinos traveling 50 Mpc to earth, assuming 
a flat distribution of relic neutrinos with a density of 54 cm -3 , is approxi- 
mately 10 28 cm~ 2 , hence the column density is not evaluated below the point 
where the cosmic rays produced on the most dilute reasonable background 
exceeds the observed flux above the GZK cut off. With the clustering of relic 
neutrinos on scales described in Ref. (0) the column density is increased 
to S(D) ~ 10 29 cm -2 , while for a neutrino cloud with radius R = 50 pc 
and density n(x) = 10 14 cm -3 we have S(D) ~ 10 34 cm -2 . This makes the 
demands on the incident neutrino flux much more modest. 

Using the fact that the flux of UHECRs is 10~ 20 < F ph < 10" 19 we can 
establish an allowed neutrino cloud parameter range, depending of course 
on the incident neutrino flux. In Sect. |2.2| . we saw that the magnitude of 
the cosmic neutrino flux is open to question. Most models place an upper 
bound of ErF(Eji,0) < 10 -19 cm -2 s~ 1 sr -1 , while the lowest expected flux is 
not at all known. However in the context of the Z-burst model we can define 
a minimum required flux, that is, the smallest neutrino flux able to produce 

1 In more familiar units ~7-7x 10 6 au or ~ 3 x — 30 pc 
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Fp/ 7 assuming the incident neutrinos are absorbed in their entirety. 



ErF v min (E R , x) = f p/7 . (7) 

ML z 

This situation would occur if the average interaction length of the neutrino 
becomes significantly less than the cloud radius. Depending on the value 
of Fp/j, 10 _20 cm _2 s _1 sr _1 or 10~ 19 cm _2 s _1 sr _1 the corresponding minimum 
neutrino flux is 5.6 x 10~ 21 or 5.6 x 10~ 20 cm _2 s _1 sr _1 respectively. The 
minimum neutrino flux is approached for high column density and small 
neutrino flux in Fig. 1. 

The neutrino cloud radius and density permitted by the observed UHE 
CR data are plotted in Figs. |2| and [3[ The neutrino flux in Fig. ^ is that 
of the upper bound E R F u (E R ,0) = 10 _19 cm _2 s _1 sr _1 , and for Fig. |3], just 
above the minimum required, ErF u (Er,0) = 10 _20 cm _2 s _1 sr _1 . In the first 
case the observed UHECRs can be fully accounted for, while in the second a 
maximum flux of F p / y ~ 1.7 x 10 _20 cm _2 s^ 1 sr~ 1 is produced. 

We now return to the question of the Fermi motion of the neutrinos in 
the cloud. This motion will alter the value of the resonant energy Er. For a 
cosmic neutrino interacting with a relic neutrino s = 2(E U E C — p v ■ p c ) where 
E c and p c are the energy and momentum of the relic neutrino, at resonance 

M 2 

and for a relic neutrino at rest 2m v E v = M§, hence E R = I n a neutrino 



cloud of high density E R will take on a range of values from E R = 

M 2 

to E R = -^-(Ef — pjcos9)~ , where Ef and pj are the Fermi energy and 
momentum of the cloud. To account for the smearing of resonant energy we 
define an average resonant energy E R , 

- 1 f p f - o Ml 



2N J ^(E c -p c cos6) 



where we have normalized to the Fermi momentum of the cloud, N = -f-Pj- 
The resulting average energy is 



E R = 4tt- 



[2p) + ml) sinh" 1 - p f jpj+ 



8N 

which in the non-relativistic limit, as expected, reduces to 



m 2 



(9) 



Mlpl Ml 

B « = ^ = ^,- < 10 > 
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In the case of a very dense cloud, with a Fermi momentum of pf ~ 13 eV/c 
and for leV/c 2 mass neutrinos, the resultant average resonant energy is E R = 
1.3 x 10 21 eV. Since our calculations have been performed in the center of 
momentum frame the reduction in resonant energy will not significantly alter 
any of our conclusions regarding the results discussed so far, in as much as the 
narrow resonance approximation holds over the width of the average energy 
distribution. 

As an aside we note that the lowering of the resonant energy also allows 
the Z-burst mechanism to remain viable should neutrinos possess very small 
masses. In a non-clustered background the resonant energy for a 0.1 eV/c 2 
mass neutrino, (rather than a 1 eV/c 2 neutrino) is E R = 4 x 10 22 eV, while 
in a dense neutrino cloud the average resonant energy can be as low as E R = 
2.3 x 10 21 eV making the Z-burst mechanism viable even if the differential 
neutrino flux drops off very rapidly with energy. 

The standard deviation of the resonant energy, both in absolute terms 
and relative to the mean energy, increases with increasing Fermi energy. 

(A£) 2 = {E R )l - {E R f. (11) 

However, for our present estimates we simply work with the mean resonant 
energy, and ignore the variation in flux over the width of the distribution. 

A stronger constraint on the parameter range of the neutrino cloud comes 
from the non-observation of coincident CR events at ultra-high energies. The 
decay products of the Z-resonance are boosted into a cone of angle 1CT 11 , thus 
occupying an opening area of A p = 7r(tan[10 _11 ]-R c ) 2 , where R c is the radius 
at which the neutrino annihilation took place. If the detection area of an 
experiment (A^ et ) is larger than A p /M then, on average a coincident event 
would be observed. The multiplicity in this context should count all particles 
produced at resonance not just those with an energy above the GZK cutoff. 

With a multiplicity of ~ 60 and a detection area of A^t ~ 500km 2 a 
sphere of radius R c ~ 10 18 cm can be defined within which neutrino annihi- 
lations will on average result in the observation of a coincident event. The 
flux of coincident events (F c ) will depend on the neutrino flux incident upon 



this sphere, and can be evaluated in the formalism of Sect. |23 
F c = 5F(E„,x>(D-R c )) 

6F(E U , D) - 6F(E U , D — R c 



-(< 

= MxE R F(E R ,0) J ^L(-exp[-a(s)S(D)] 

+ exp[-a(s)S(D - R c )} ) . (12) 
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If in addition we make the simplifying assumption that the relic neutrino 
column density outside the cloud is much less than the column density inside 
we are able to set n(x) = n c (x). Now the expression for the flux of coincident 
events becomes 

r ds 

F c = E R F(E R ,0) J — exp[-a(s)n c R] 

x ( exp[a(s)n c R c ] - l) . (13) 



The non-observation of a coincident flux implies the following bound on 
the total number (N t ) of such events observed by a particular experiment 

N t = A t r F c < 1, (14) 

where A is the experimental aperture at Eq > 10 19,6 eV and t r is the experi- 
mental running time. In order for the neutrino cloud hypothesis to hold true, 
values of cloud density and radius must be found so that for each experiment, 

r ds 

F c > E r F(E r ,0) / — exp[-a(s)n v R) 

x ( 1 - exp[a(s)n u R c ] ) . (15) 



The AGASA experiment currently has the largest exposure at 

At r = 670 km 2 sr year ~2x 10 21 cm 2 sr s . (16) 

This allows a coincident flux of F c ~ 5 x 10 _22 cm~ 2 s _1 sr _1 . The values of 
R and n v which solve Eq. [15| are shown in Figs. @ and || The coincidence 
constraint allows a minimum cloud radius of -R m i n = 10 18 cm to be specified, 
and in some sense a maximum relic neutrino density, depending on the inci- 
dent neutrino flux. The results of this calculation, and the bounds from the 
UHECR data serve to define an allowed parameter region for the neutrino 
cloud; see the shaded area in Figs. 2 and 3. We emphasize that, while the 
general principle that the non-observation of coincident events will eliminate 
a range of cloud radii less that a some "minimum" radius, the particular 
value of this minimum radius is determined only after a model for the cloud, 
is chosen, and the number density is calculated. The uniform cloud consid- 
ered here is but the simplest of the cloud models introduced in Ref. (0), 
and even in that paper only a few simple cases were considered. These other 
models, with considerations similar to the above, will lead to different values 
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for the minimum radius of the cloud. We do not pursue this discussion into 
limitations regarding the variety of possible clouds in the present paper. 

The constraints on the parameter range of neutrino clouds as a result of 
the arguments based on coincident events have not taken into account the 
spatial resolution of the various CR experiments. While this is not an issue 
for medium to large clouds, it does become important for clouds with smaller 
radii. For example, if a neutrino cloud had a radius of 50au decay products 
would be confined to a cone of radius less than 10m at earth. This radius 
would be too small for most CR experiments to resolve, in which case the 
decay products will be measured as a single event with an energy of E R . 
In the minimal Z-burst model Er is a significantly higher energy than the 
highest energy CR event. However as we have shown in this note, for high 
density clouds the effect of the Fermi motion is to reduce the average value 
of E R so that the possibility of small radii clouds is not excluded. We leave a 
full quantitative analysis to future work, and for the time being acknowledge 
that there is a region of allowed parameter space for small radii neutrino 
clouds. 

4 Conclusion 

We have shown it is possible that neutrino clouds, interacting with UHE 
neutrinos, are capable of producing the measured flux of UHECR's, and that 
this result is within both theoretical and experimental bounds imposed on 
the incident neutrino flux. The dimensions of the neutrino clouds required 
to produce the measured UHE CR flux assuming the Z-burst mechanism is 
responsible were investigated, and need no new theoretical treatment beyond 
that given in Ref. (0). The requirement that no coincident events are ob- 
served requires that uniform clouds, larger than some small radius still to 
be determined, are at least 0.3pc in radius — rather larger than originally 
envisaged in Ref. (0), but without disturbing the viability of the concept 
of neutrino clouds. Clouds with different density profiles will have different 
limitations placed on their radii. 
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Figure 1: The dashed line corresponds to the required column density needed 
to produce a photon/hadron flux of 10 _20 cm _2 s _1 sr _1 , while the solid line is 
the column density need for a photon/hadron flux of 10 _19 cm _2 s _1 sr _1 . The 
shaded region corresponds to the experimentally allowed column density. 
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Figure 2: The solid line represents the cloud parameters required to produce 
a GZK flux of -F p / 7 = 10~ 20 cm _2 s _1 sr _1 , with an incident neutrino flux of 
ErF u (Er, 0) = 10 _19 cm~ 2 s _1 sr _1 . The dashed line represents a GZK flux of 
F p / 7 = 10 _19 cm _2 s _1 sr _1 . The short dashed line is the bound on R result- 
ing from the non-observation of coincident events with an incident neutrino 
flux of E R F v (E R ,ti) = 10 _19 cm~ 2 s _1 sr _1 . The shaded region is the allowed 
parameter space. 
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Figure 3: The solid line represents the cloud parameters required to produce 
a GZK flux of F p / 7 = 10~ 20 cm _2 s _1 sr -1 , with an incident neutrino flux of 
E R F u (E R ,0) = lO-^cm-V 1 ^- 1 . The dashed line represents a GZK flux 
of .F p / 7 = 1.7 x 10 _20 cm _2 s _1 sr _1 , the maximum flux for this small inci- 
dent neutrino flux. The short dashed line is the bound on R resulting from 
the non-observation of coincident events with an incident neutrino flux of 
E R F y (E R) Q) = 10 _20 cm _2 s _1 sr _1 . The shaded region is the allowed param- 
eter space. 
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